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ANALYTICAL CHEMICAL KINETIC INVESTIGATION OF THE EFFECTS 
OF OXYGEN, HYDROGEN, AND HYDROXYL RADICALS 
ON HYDROGEN -AIR COMBUSTION 

By George T. Carson, Jr. 

Langley Research Center 

SUMMARY 

Quantitative values have been computed which show the effects of oxygen, hydrogen, 
and hydroxyl radicals on the finite-rate chemical kinetics of premixed hydrogen and air 
undergoing isobaric combustion. Mole fractions equal to 10 10 "5, and 10 of the 
free radicals were considered to be initially present, either singly or in combination, in 
homogeneous, hydrogen -air mixtures with equivalence ratios of 0.2, 0.6, 1.0, and 1.2. 
Initial temperatures were 1100 K, 1200 K, and 1500 K and pressures were 0.5, 1.0, 2.0, 
and 4.0 atm. A few other initial states were used when necessary to determine inter- 
mediate trends. The investigation demonstrates that small amounts of radicals reduce 
induction time, herein defined as the time to 5 percent of the temperature rise from the 
initial mixture temperature to the calculated equilibrium temperature. Also, the study 
shows that within the range of conditions investigated, hydroperoxyl (HO2) plays a sig- 
nificant role in the early part of the ignition process; the inflection in its rate of forma- 
tion can be used to define induction time. The reaction time, defined as time between 
5 percent and 95 percent of the temperature rise, is not affected by the initial amounts of 
free radicals. All of the radicals are more effective in reducing induction time at low 
initial temperatures. For all mixtures and states, atomic oxygen is the most effective 
radical of the radicals investigated. The results represent equilibrium chemistry, 
finite -rate chemical kinetics, and curve -fitting computer programs over this range of 
conditions. It is also implied that additives which yield free radicals can be used to 
effect a compact combustor design or to permit operation in otherwise reaction -limited 
regimes of initial pressures and temperatures. 

INTRODUCTION 

In support of the NASA Hypersonic Research Engine Project, discussed in refer- 
ence I, an analytical chemical kinetic investigation of the effects of free radicals on 
hydrogen-air combustion was conducted. The most difficult problem areas in the design 
of the hypersonic ramjet engines are associated with the combustion process. The flow 



through the combustor Is continuous and supersonic and is composed of mixing and of 
burning fuel and air. The fuel injection and mixing parts of the processes are extremely 
important, and in diffusion controlled combustion they control heat release. Where the 
mixture of fuel and air at the beginning of the combustor is near the threshold tempera- 
ture for autoignition, the combustion kinetics plays a major part in ignition. In order 
to minimize the combustor length (thus lowering weight and overall heat transfer), 
shortening the time for completion of the combustion reaction is critically important. 

The initial testing of supersonic combustion ramjet engines will be in ground-based 
facilities. The high -temperature air required for simulation can be obtained by first 
heating nitrogen and then mixing it with oxygen. Another method of producing a hot test 
gas is to burn fuel in a high-pressure oxygen -rich airstream. The combustion process 
in the burner, or "heater,” produces a hot -gas flow whose composition has the proper 
percentage of oxygen to simulate air. The main difficulty with the latter method is that 
the air ingested by the test engine contains combustion products, as well as unrecom- 
bined free radicals; this may produce misleading test results. 

The time for combustion (i.e., the total time associated with completion of the 
chemistry of combustion) is usually divided into two parts: induction time and reaction 
time. For induction time, investigators sometimes use a criterion based on reaching a 
prescribed concentration level of some product species (ref. 2 uses the hydroxyl radical 
at the level of 10-6 mole per liter), or one based on the rate of temperature rise due to 
three -body radical recombination (ref. 3), or one based on experimentally measured 
light emission or absorption (ref. 4). In the present paper, as in reference 5, induction 
time is defined as the time to reach 5 percent of the equilibrium temperature rise, i.e., 

5 percent of the difference between the initial static temperature of the mixture and the 
equilibrium temperature. The reaction time is the remainder of the time required to 
achieve complete combustion. However, since the chemical process approaches com- 
pletion asymptotically, the latter part of the heat release is lost in the exhaust of a 
scramjet. Therefore, for the purpose of this paper, a finite limit for the reaction was 
chosen to be 95 percent of the potential temperature rise. Reference 5 includes a lim- 
ited comparison of the induction times of hydrogen -air mixtures with and without initial 
amounts of the free radicals, oxygen (O), hydrogen (H), and hydroxyl (OH). It is shown 
that the presence of initially small amounts of free radicals results in a significant 
reduction of induction time. The present analysis was undertaken in an attempt to 
improve the understanding of the significance of small amounts of contaminants in com- 
bustible gaseous mixtures. 

Only the chemical kinetics of a premixed one -dimensional stream of hydrogen and 
air is considered here. Elimination of injection and mixing was considered a necessary 
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simplification in order to achieve the objectives. The chemical kinetic analyses of 
hydrogen-air combustion at constant pressure (diverging combustor) which are presented 
show the effects of small amounts of atomic oxygen, atomic hydrogen, or hydroxyl radi- 
cal on the Induction and reaction times. The range of thermodynamic state points and 
compositions included are typical of ramjet operation. The scope of the computer pro- 
grams used in this analysis precluded the investigation of mixing length and local equiv- 
alence ratio variance. 
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SYMBOLS 

number of chemical species 

specific chemical species 

finite -rate reaction 

reaction rate constant for reaction j 

reaction rate constant for reverse of reaction j 

molecular mass of species i 

static pressure (1 atm = 101 325 Pa) 

universal gas constant, 8.31434 J/mol-K 

static temperature, K 

equilibrium temperature, K 

mean static temperature, ^(Tf + Tf), K 

initial static temperature, K 
time, psec 

effective mole fraction of radical species 
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Xi 

Yi 

At 

At/tq 

9 

T 

'^O 


mole fraction of specific chemical species, 


Moles of species i 
Moles of mixture 


mass fraction of specific chemical species, 


Mass of species i 
Mass of mixture 


induction -time difference, Tq - r, psec 

relative reduction in induction time 

fractional temperature rise, (t - " "^l) 

induction time of a mixture (time at which 9 = 0.05), fisec 

induction time of a mixture not containing initial quantities of free radicals 
O, H, or OH, psec 

reaction time (time to 6 = 0.95 less time to 0 = 0.05), psec 
react ion -time parameter, psec-(atm)^‘®^^ 


0 


equivalence ratio, 


Fuel -air ratio 

Stoichiometric fuel -air ratio 


ANALYSIS 


In order to obtain numerical values used in this analysis, two different computer 
programs were used. First, the equilibrium chemistry program (ref. 5) determines the 
equilibrium temperature and composition for mixtures at given initial temperatures and 
composition under isobaric conditions. Second, the chemical kinetics program (also by 
the authors of ref. 5) determines temperature and composition in mixtures reacting at 
finite rates. An isobaric, premixed, one -dimensional adiabatic flow is assumed. 

Four classes of gaseous mixtures were considered: 

(1) Air: As a basic simplifying assumption, assume air to consist only of nitrogen 
and oxygen; this seems justified since other gases are present only in small amounts and 
the radicals to be investigated are not present. Reference 6 indicates that air contains 
20.9472 -percent oxygen and 79.0528 -percent nitrogen by volume. Also, since the molec- 
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ular weight of air is given as a constant to an altitude of 90 km in reference 6, the com- 
position of air is assumed to be constant. 

(2) Stoichiometric fuel-air mixture: A stoichiometric mixture of fuel and air is 
one in which both reactants can react completely to form products. The stoichiometric 
hydrogen -oxygen reaction is indicated by the equation 

2H2 + 02 = 2 H 2 O 


The oxygen -hydrogen ratio in a stoichiometric mixture is determined by the relationship 

. XiMi 
Y, --J- 

which gives 

YOo 

7.936507 

The nitrogen-oxygen mass -fraction ratio for air is similarly determined as 

-^= 3.3103448 
^02 


From the definition of mass fraction it is known that 

I 

^ Yi = YN2.Yo2 + YH2=l 
i=l 


The three unknown quantities can be found from the preceding three equations. The 
formation of oxides of nitrogen below 3000 K for short residence times is assumed 
negligible. 

(3) Nonstoichiometric fuel -air mixtures: The mass fractions of nitrogen, oxygen, 
and hydrogen for nonstoichiometric compositions are computed in a similar manner. 

The difference is the inclusion of an equivalence ratio <j) which in this report is defined 


as 


0 = 


C^H2/^02)stoichi 


stoichiometric 
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An equivalence ratio of unity signifies stoichiometric proportions, less than unity implies 
a fuel -lean mixture, and greater than unity implies a fuel -rich mixture. The composi- 
tion of baseline mixtures in mass fractions is shown in table 1. 

(4) Small amounts of initial free radicals: Finally, both stoichiometric and non- 
stoichiometric mixtures were calculated that had small (typically, 10"^ to 10"^) mole 
fractions of free radicals, O, H, and OH, initially present. 

The chemical reactions and rate constants kj in this analysis are shown in 
table 2, along with their sources (refs. 7 to 16). Although rate constants and third-body 
efficiencies used in the chemical kinetics program may be easily changed, the reactions 
themselves cannot be changed without reprograming. For the investigation this was not 
considered to be a drawback because the relative effect of free radicals was of prime 
importance. Since the effects of the free -radical reactions with hydrogen peroxide and 
hydroperoxyl are small above 1100 K, they were not included in the analysis. At the 
time of the preparation of the present paper, the authors of reference 16 experimentally 
determined a new value of k]^2> value is designated herein as k^2A* 

third-body efficiencies are given in reference 5. The value for reaction 8 is set at unity 
in order to establish a baseline reacting mixture which does not depend upon water as a 
third body. 

The first parameter derived from the computational values was a fractional tem- 
perature rise which was defined as 


Tf - Ti 

This parameter expresses quantitatively the amount of reaction corresponding to a 
change in the mixture temperature. In this investigation, as in reference 5, induction 
time T is defined as the time at which 9 = 0.05 and the total reaction time is defined 
by 6 s= 0.95. Another parameter used is the relative reduction of induction time. It 
has been given the symbol At/tq and is defined as the ratio of the difference between 
the induction time of a mixture with and without initial free radicals At = r - Tq to the 
induction time without initial free radicals. 

The computed values of the relative reduction of induction time and their relation- 
ships with the prime parameters of this study are discussed in the following section. 

RESULTS AND DISCUSSION 
Equilibrium Chemistry Computations 

The results of the computations, as shown in figures 1 to 4, give equilibrium tem- 
perature as a function of initial temperature along isobars for each of the equivalence 
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ratios. These curves are valid with or without initial amounts of free radicals up to 
10 mole per mole of mixture. 

Finite -Rate Chemistry Computations 

Temperature -time mapping. - Figures 5 to 8 map static temperature as a function 
of time for a constant-pressure combustion process. Shown are lines of constant initial 
temperature and fractional temperature rise. These figures demonstrate the impor- 
tance of initial temperature. Below 1000 K the induction time becomes so great that the 
fuel -air mixture cannot ignite during the residence time of about 250 fi sec. Figures 5, 

6, 7, and 8 represent examples at which the static pressures p are 0.5, 1.0, 2.0, 
and 4.0 atm, respectively; and comparison of the data of these figures shows that 
increasing pressure decreases induction time. This seems to be generally true except 
for low initial static temperature (i.e., near 1100 K) where increasing pressure has an 
inhibiting effect, as is clearly depicted by figures 7 and 8. For T]^ = 1100 K, 
p = 4.0 atm, <p = 1.0, and 6 = 0.5, ignition did not occur even beyond 1000 psec. 

Table 3 is a presentation of input values T, p, Xq, ^Oh) 

resulting induction times and related fractions (r, t/tq, and At/tq). Since atomic 
oxj^en produced the most dramatic results, it was used for graphic presentation. 

Figure 9 shows the effect a 10"^ mole fraction of atomic oxygen has on combustion. 
In comparing this figure with figure 6, it is noted that the presence of initial radicals 
has an acceleration effect on combustion which is large for low initial temperature but 
is negligible for high initial temperature. Because of the large range in magnitude of 
the parameter, graphs subsequent to figure 9 are plotted on logarithmic scales. 

Induction -time mapping.- Figures 10, 11, and 12 show the variation of induction 
time with static pressure for families of initial temperature, equivalence ratio, and ini- 
tial trace amount of radical, respectively. From figure 10 it is apparent that the effect 
of initial temperature is large. For initial temperatures of 1200 K (and below) there are 
definite minimum induction times. For 1500 K and above, no minimum occurs below 
10 atm. According to reference 15 this is because reaction 9 inhibits the hydrogen -air 
reaction at low temperature and high pressure. The effect of equivalence ratio is minor, 
as may be observed from figure 11. The pronounced effect of small amounts (10 10 "5^ 
and 10-^) of initial atomic oxygen is shown in figure 12. 

^ Time variation of species . - The variation of species mass fraction with time is 

shown in figures 13 to 18. Figure 13 depicts a stoichiometric mixture of hydrogen and 
air initially at 1100 K and without any initial radicals. The increase of products (radi- 
cals and H 2 O) is very rapid. Figure 14 shows the results of adding small amounts of 
radicals and indicates that the formation of products begins much sooner when the radi- 
cals are present. However, the rate of change of species mass fraction with time is 
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much less. The chemical -kinetics prc^ram shows that the reactions containing atomic 
oxygen begin at a low temperature and, being initially of small magnitude, produce little 
temperature rise to increase the reaction rates. Comparison of figures 13, 15, and 17 
shows that, for increasing initial temperature (i.e., 1100, 1200, and 1500 K, respectively), 
the products form sooner even though the rate of change is less; also, when the compar- 
ison is expanded to include figures 14, 16, and 18, the product -forming action of atomic 
oxygen becomes less pronounced as the initial temperature is increased. In fig- 
ures 14(a), 16(a), and 18(a), which cover the first part of the reaction only, it is observed 
that the amount of atomic oxygen actually decreases at first. This pheonmenon is appar- 
ently because of reaction 7 which results in removal of atomic oxygen and hydrogen. 
Figures 13(a), 15(a), and 17(a) show that all products monotonically increase for mix- 
tures which do not contain initial atomic oxygen. 

Hydroperoxyl relation to induction time , - It was observed that the peak of hydro - 
peroxyl formation occurred near the selected Induction time criterion of 0 = 0.05. This 
trend was investigated for cases which have large induction times, i.e., low initial tem- 
perature and high pressure. Figures 19 and 20 show that this trend holds. Therefore, 
it appears that the peak of hydroperoxyl formation is a good indicator for finding the 
induction time. For the range of temperature and pressure investigated, the two equa- 
tions shown in table 2 involving hydroperoxyl appear to be sufficient. 

Figure 21 shows the relationship between induction time and static pressure as 
calculated with and without the formation of hydroperoxyl. It is evident from the figure 
that computations which do not include hydroperoxyl give increasingly erroneous results 
as the static pressure increases. According to reference 17, HO2 is formed in chain- 
breaking reactions, i.e,, reactions which remove radicals from the mixture, thereby 
inhibiting product formation. 

Relative reduction in induction time . - Figures 22 to 25 illustrate the relative 
reduction in induction time as a function of mole fraction of initial radical at pressures 
of 0.5, 1.0, 2.0, and 4.0 atm. These figures show that the effectiveness of O, H, and OH 
radicals in reducing induction time is inversely related to initial temperature and tends 
to level out as mole fraction increases. Atomic oxygen is consistently more effective 
than atomic hydrogen and hydroxyl. Pressure has an effect on the relative reduction in 
induction time, but it is not consistent. The omission from figure 25 of a relationship 
for the initial temperature of 1100 K is due to the fact that the induction time for these 
state conditions is excessive, i.e., beyond 400 jisec, as may be seen in figure 8. There- 
fore, At/tq was not calculated. 

The peculiar relation between At/tq and pressure at various atomic oxygen mole 
fractions is presented in figure 26. Although the trends show some measure of depend- 
ence, the effect of pressure is minor below 2.0 atm. At 4.0 atm the effect is pronounced. 
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For values of Xq larger than 0.3 x 10“^, pressure reduces induction time; for Xq 
less than 0.3 x 10~^, pressure increases induction time. This phenomenon is presumed 
to be caused by third-body reactions which involve atomic oxygen. When the initial 
amount of this radical is present in sufficent quantity, it interacts with other species and 
thereby enhances reaction. However, when the quantity of atomic oxygen is very small, 
the pressure dependent reaction 6 becomes predominant and inhibits the reaction. 

The next parametric investigation was the effect of (p on At/tq for a family 
of values of Xq. The results are displayed at two pressures, 1.0 and 4.0 atm, in fig- 
ures 27 and 28, respectively. The variation of At/tq with (p is not particularly 
remarkable; this is as expected, since even at stoichiometric proportions the mass frac- 
tion of H 2 is much less than for N 2 or O 2 . In figures 27 and 28 it is noted that the curve 
for Xo = 10"^ does show an inverse relationship of At/tq with (p. 

The variation of At/Tq with Tj for a family of Xq is depicted by figure 29. 

As is roughly indicated in figures 5 to 9, At/Tq is inversely proportional to Tj, also, 
from figure 12 and figures 22 to 26, it can be seen that increasing values of Xq cor- 
respond to increasing values of At/Tq. An attempt was made to determine a parametric 
expression from the information generated by the computations. The plan was to find 
coefficients of Xq, Xjj, and Xqjj which would correlate the curves of At/tq ver- 
sus radical species mole fraction into a single line for each Tj^. Next would be corre- 
lation of initial temperatures by finding the required function of Tj^. Finally, the effects 
of p and (p would be included. 

An abscissa transformation performed for the curves for Xjj and Xqjj relative 
to the curve for Xq on each of figures 22 to 25; it was determined that H and OH were 
approximately 60 percent as effective as O in increasing At/tq. Therefore, an effec- 
tive mole fraction was defined as 

Xe = Xq + 0.6 (Xjj + Xqjj) 

Since the 60 -percent approximate effectiveness appeared valid for 0.5, 1.0, and 2.0 atm, 
the effects of pressure were considered masked. It should be noted that the 60 -percent 
value applies to the set of rate constants used. 

Merging of the initial temperature family of curves into a single expression was 
not successful, although curves expressing At/tq versus Xg were found for the 
intermediate values of Tj = 1300 K and 1400 K by linear interpolation between 
Tj = 1200 K and 1500 K. Figure 30 shows a family of Tj curves, including the inter- 
polated ones of 1300 K and 1400 K. 

It is noted that the lack of accuracy of rate constants and third-body efficiencies 
form a limit to the accuracy of this analysis. As new experimental values are deter - 
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mined, they may easily be substituted into the program. For example, the rate constants 
for reactions 2 and 3 were changed and are listed in table 2 as k 2 ^ and Consi- 

der two mixtures at an initial temperature of 1200 K, a pressure of 1 atm, and an equiv- 
alence ratio of unity, one mixture without initial free radicals and one mixture with an 
initial mole fraction of 10 of atomic oxygen. The rate constant k 2 A gave induction 
times of 38 Msec and 14,5 Msec for the two mixtures, respectively, and the rate constant 
^3 A induction times of 41.86 Msec and 16.1 Msec for the two mixtures, respec- 

tively, compared with the values of 40.4 Msec and 15.32 Msec obtained using k 2 and 
k 3 . Also, the relative reductions in induction time for k2^ and k 3 ^ were found to 
be 0.618 and 0.615, respectively, compared with 0.621 using k 2 and k 3 . 


Before this investigation was completed, an updated experimental value for the 
rate of the twelfth reaction (H 2 + O 2 = OH + Oh) became available from reference 16. 

The newer one (designated 12A) was used in the computation of induction time. Figure 31 
is a comparison of induction time relative to initial temperature for computations using 
ki2 and kj2A- ^12A ^e^^^ces induction time, but only very slightly. 

In order to compare the criteria for the induction time of references 2 and 5, the 
concentration of hydroxyl radical was calculated for two mixtures. Both mixtures had 
an initial temperature of 1200 K, pressure of 1 atm, and an equivalence ratio of unity. 
However, one mixture initially contained 10“^ mole fraction of atomic oxygen. When the 
mixtures had experienced 5 percent of the equilibrium temperature rise (ref. 5), the 
hydroxyl concentration was found to be 1.4 x 10“^ mole per liter and 1.8 x 10"^ mole per 
liter for the mixtures without and with initial atomic oxygen, respectively. This was 
believed to be a reasonable agreement with the criterion of reference 2. Also, the cri- 
teria for references 3 and 5 were compared. Using the same mixtures as before, it was 
found that maximum atomic hydrogen concentration, i.e., before beginning of third -body 
recombination (ref. 3), occurred at 48 Msec and 23 Msec for the mixtures without and 
with atomic oxygen, respectively. This compared fairly well with 40.4 Msec and 
15.32 Msec for the induction time as defined by the criterion of reference 5 and used in 
the present paper. 


Reaction time . - The final portion of this investigation considered the effects of 
trace amounts of O, H, and OH radicals on reaction time. The reaction time was found 
to be unaffected. Figure 32 is a map of reaction time as a function of the reciprocal of 
mean temperature for intersecting families of pressure and Initial temperature. The 
mean temperature was used in order to obtain a direct comparison to reference 5 as 
shown in the figure. Reference 5 indicated a correlation given by 


log = M 8 0 + 0.040 

Tm 
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Several of the rate constants are different from those used in reference 5, and the 
third-body efficiency of reaction 8 is four times as great. The value of third -body effi- 
ciency did not affect the induction time but was inversely related to the total time, 
thereby changing the reaction time. Thus, the accuracy of the analysis is limited by the 
accuracy of the rate constants and the third-body efficiencies. The missing regions of. 
the graph represent those low initial temperature, high pressure regimes for which the 
induction times were too large to be calculated. The parallelism of the isobars on the 
small semilogarithmic plot si^gested that they could be condensed to a single curve. 

The curve obtained by correlation was of the form 

log + 0.47 

Tm 

which plots coincident with the curve p = 1.0 atm and the ordinate expressed as 
lasec-(atm)^-^^^. 

A more convenient independent variable for reaction time is initial temperature. 
This relationship is shown by figure 33. Also shown in this figure is a comparison from 
reference 18. The correlation given by reference 18 is 

^ 1000 

which can also be expressed as 

log = -0.0004486T1 + 2.021 

This correlation determined in the present analysis was 

log Tj.pl‘641 = -0.000415Ti + 2.427 

CONCLUDING REMARKS 


Quantitative analytical results have been computed which show that small amounts 
of atomic oxygen, atomic hydrogen, and hydroxyl radicals reduce induction time of 
hydrogen -air mixtures undergoing autoignition. 

The peak of the reaction chain-breaking species, hydroperoxyl, is especially defin- 
itive and seems to be a logical criterion for defining induction time. For mixtures either 
containing or not containing initial free radicals, the peak of hydroperoxyl occurs con- 
sistently near the time of 5 -percent fractional temperature rise. 

' Relative reduction in induction time is found to be a convenient index to measure 
the comparative effectiveness of the different initial free radicals, as well as their con- 
centration, Atomic oxygen was consistently more effective than atomic hydrogen or 
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hydroxyl. On the average, both atomic hydrogen and hydroxyl were found to be 60 per- 
cent as effective as atomic oxygen. 

The equilibrium temperature was not affected by the presence of small amounts of 
radicals. Evidently, this was due to the fact that the initial and equilibrium composition 
and enthalpy values were not significantly changed. 

The addition of small amounts of radicals to the initial reacting gases did not 
affect the reaction time (the time between 5 and 95 percent of the equilibrium 
temperature). 

The accuracy of this analysis was limited by the accuracy of available reaction- 
rate constants and third -body efficiencies. Both of these quantities are easily changed 
in the chemical kinetics computer program that was used. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 19, 1974. 
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TABLE 1.- COMPOSITION OF BASELINE MIXTURES 
IN MASS FRACTIONS 


<p 

Yn2 

O 

to 


0 

0.768 

0.232 

0 

.2 

.7635343 

.230651 

.0058124 

.6 

.7547621 

.2280011 

.0172367 

1.0 

.7462104 

.2254107 

.0284017 

1.2 

.7419728 

.2241376 

,0338896 



TABLE 2.- REACTION HATE CONSTANTS 


j 

Reaction 

Reaction rate constant 
(a) 

Source 

1 

OH + Hg g— ^ H 2 O + H 

ki = 2.3 K 10 l 3 exp(-520Q/RT} 

Dixon -Lewis, Wilson, and Westenberg (ref. 7} 

2 

H + 0 ? OH + 0 
k _2 

k 2 = 9.9 X 10^^ exp(-15 020/RT) 

Jachimowskl and Houghton (ref. 8 ) 

2A 

H + O 2 ^ OH + 0 
*^- 2 A 

k 2 A = 2.24 X exp{-16 800/RT) 

(ref. 9) 

3 

0 + Hn izS OH + H 

kg = 7.5 X 10^3 exp(-ll 100/RT) 

Jachlmowski and Houghton (ref. 8 ) 

3A 

^3A 

0 + H? OH + H 

kjA = 4.3 X 10^3 exp(-10 200/RT) 

Wong and Potter (ref. 10) 

4 

OH + OH T — H 9 O + 0 

k^ = 1.5G X 10 I 2 exp(1000/RT) 

Dixon -Lewis, Wilson, and Westenberg (ref, 7) 

5 

^5 

H+H + mZZ 1 h 2 +M 
^-5 

kg 8.5632 X 1 d 11 t 0-45 exp(U 590/RT} 

Myerson, Watt, and Joseph (ref. 11) 

6 

kg 

0 + 0 + M O 2 + M 

k ,6 

kg 02 = 1 d 1 %- 1'22 

Kiefer and Lutz (ref. 12) 

7 

krj 

O+H + M^OH + M 

ky^lVI = ^‘0 

Kaskan and Browne (ref. 13) 

8 

H + OH + M HgO + M 

= 4.5 X 1o21t-1-5 

Kaskan and Browne (ref, 13) 

9 

H + O 2 + M ;:r HO 2 + M 

kg = 2.0 X 10^^ exp(870/RT) 

Gutman, Hardwidge, Doughert, and Lutz (ref. 14) 

10 

^10 

HO 2 + Hg ^ H 2 O 2 + H 

k -10 

kjQ = 9.6 X I 0 I 2 exp(-24 000/RT) 

(ref. 9) 

11 

^11 

H 2 O 2 + M OH + OH 4 M 

•^11, Ng = ^ exp(-46 300/RT) 

Baldwin and Bratton (ref. 15) 

12 

^12 

H 2 + O 2 ^ OH 4 OH 
^-12 

ki 2 = 1.05 X 10^^ exp(-48 080/RT) 

Jachitnowski and Houghton (preliminary to ref. 16) 

12A 

^12A 

H 2 + O 2 OH 4 OH 

’^-12A 

kj2A = exp(-46 150/RT) 

Jachimowskl and Houghton (ref. 16) 


^All reactions have units of cmS/mol-sec except reactions 5 to 8, 10, and 12 which have units of cm^/mol^-sec; activation energies 
are expressed in units of cai/mol (1 cal/mol = 4.184 J/mol). 
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Figure 4.- Equilibrium temperature as a function of initial temperature along isobars for (p = 1.2. 
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Figure 7.- Isobaric combustion temperature variation with time for p = 2.0 atm and <p ~ 1.0. 
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Figure 9.- Isobaric combustion temperature variation with time for p = 1.0 atm, (f> = 1.0, and 
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Figure 10.- Induction -time variation with static pressure for a family of initial temperatures 

for 0 = 1.0 and Xq = 0. 
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Figure 11. > Induction -time variation with static pressure for a family of equivalence ratios 

for Ti= 1200 K and Xq=0. 
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Figure 12.- Induction -time variation with static pressure for a family of initial mole fractions of atomic oxygen 

for Tj = 1200 K and 0 = 1.0. 
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(a) Mass fraction, 10"'^ to 10“'^. 


Figure 13.- Time variation of species mass fraction for = 1100 K, p = 1.0 atm, 

4 > - 1.0, and Xq = 0. 
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Figure 13.- Concluded, 
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Figure 14.- Time variation of species mass fraction for = 1100 K, p - 1.0 atm, (p - 1.0, and Xq:= 10"^. 


Species mass fraction 


Equilibrium 



Time, t, ^sec 

(b) Mass fraction, lO""^ to 1. 
Figure 14.- Concluded. 
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(a) Mass fraction, 10 to 10“^. 

Figure 15.- Time variation of species mass fraction for T j = 1200 K, p = 1,0 atm, 

4 > = 1.0, and Xq = 0. 
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Figure 15.- Concluded. 
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Figure 16.- Time variation of species mass fraction for Tj = 1200 K, p = 1.0 atm 
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Figure 16.- Concluded. 


44 


Species mass fraction, 


iMBaaBBaiiiKKsaBaEaaTiaaNaaHaBBBetfisaiBHKimafiaamBBieagBBgKBagraBaaaaagMiiiiHBnBMWHUHHiMeaBagil 


■BgnBgMiim?BaiBimK:BiHBMaHHmgaHmBaMaiBH.aiB«BiHiiBaaHiga3HHa 


s:B;:::&EU8:sasifflr!i!9B;flsaE&sB:BHi«aG:aB‘.9BinRBisnBSRaB5sa8HBiiinB!8us:ss:»ii:gi::s:s 


.nBUBOIUIOUH»rO»jrjnH»ff?IBrJlBfflfflilBL«I^BSnilll3:SIDBBH!SSHfflKI 


irissii 


10' 



SiiiSaBi^Ei^iii 


IiS!| 


MiB ttBiBSflPj 

mnEBifli 


straEssaauBj 
MJKaiRBESB&SS; 


l*j«SiH!!SEsB5ilKl4M 
lilEicaauBBasai 

lilSB 


Iliisii 


aaBBEBi 

aSIBBSSSB 


-'/MraBnBnHnnMi 

HFiiSBfflBiaiii&ssHn ass 


orsaiomaiissaseiai 

liSBESSfiB 




liiiii 

g'B3QmBaaiaaiu:sBs:i’«BaBBaKB8S»BiiHiffiisBBUKnBUR]BB:aui: 

VRBSimiBaiEaflHrHHBBBnHRiRailllRBBBSaiinUIIUBBBn 

^MiBHgHnaiHimggniagB-gggggaKHmaBBiHgiiniBiaBBBaK 


:IR»::!HSSS8EBK 


iniia 


9SS! niilsHlaaiSEiu 
^H!»BBSBB8E: 


JiiirBRuniiiaEniiiHrMiifiBffiinHBnHB 

rirsRgesr^iiiiiieii^^^Beg^i 


SRiSU 


aSniBS 

snuum 


uuiiBauHgil 


li Biiijia i 

U ■■■■11 I I 


igaHm 

»!BI 


I II ■■■s 

:*t aaeBBB 

iEumfliEil 


!;»|BSaS »KSBBBSB»SE!UKIBU8BBI 


ipisai 

siS 

ISU, 

ail HI 
ail mi 



Eiilhlliili 


iBBSSiiigiiHslii 


inmwBmillBSIHinmUBBSinUBBBnKBSISiailBI&SSSEanilBUIBBBBBSR 

siBiaaKiRnfss:»RinisnraBiiBnE»ssssfl:::E%^Biss:s:siaii3flSBS3SR, 

l8BUSS{R?n’l»:9ilBm'.ia6L^K8Sn8BHe58SS:gHll!£l'.aHni6MlBUIIII 

iBHSUKii8S!snu;uiPKsiss»«HBB:ss::BiURi&TisHaRR}«iii9ai:nBi 

lilMi^eigg!SiiiiS5T=iillsr 


llpigsili 

iniiBiasa 


ISICSEmBB 

«S=S83B 

mm 


4 6 8 


4 6 8 jq2 


4 6 8 


Time, t, p,sec 

(a) Mass fraction, 10“"^ to 10"'^. 

Figure 17.- Time variation of species mass fraction for Tj = 1500 K, p = 1.0 atm, 

<f> - 1,0, and Xq = 0. 
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Figure 17.- Concluded. 
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Species mass fraction, 



Figure 18.- Time variation of species mass fraction for Tj = 1500 K, p = 1.0 atm, <p = 1-0, and Xq = 10~4. 
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(b) Mass fraction, 10"'^ to 1. 
Figure 18.- Concluded. 
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Species mass fraction. 



Figure 19.- Time variation of the mass fraction of H 2 O, HO 2 , 
H 2 O 2 for = 1100 K, p - 2.0 atm, and <p ^ 1.0. 
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Fi^re 20.- Time variation of the mass fraction of H2O, HO2, 
and H2O for Tj = 1100 K, p = 4.O atm, and <p = 1.0. 
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Figure 21.- Induction time variation with static pressure, considering and not considering hydroperoxyl 

for Ti = 1200 K and 0 = 1.0. 
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Figure 22.- Relative reduction in induction time as a function of mole fraction of atomic 
oxygen, atomic hydrogen, and hydroxyl for p = 0.5 atm and = 1.0. 
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Relative reduction in induction time, At/t^ 
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Figure 23.- Relative reduction in induction time as a function of mole fraction of atomic oxygen, atomic 

hydrogen, and hydroxyl for p = 1.0 atm and <p ~ 1.0. 
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Figure 24.- Relative reduction in induction time as a function of mole fraction of atomic 
oxygen, atomic hydrogen, and hydroxyl for p = 2.0 atm and </> = 1.0. 
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Figure 25. - Relative reduction in induction time as a function of mole fraction of atomic 
oxygen, atomic hydrogen, and hydroxyl for p = 4.0 atm and 0 = 1.0. 



55 


Fjftlative reduction in induction time, At/t 



Pressure, p, atm 

Figure 26.- Variation of relative reduction in induction time with pressure 
for Ti = 1200 K and 0 = 1.0. 
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Relative reduction in induction time, At/t 



Figure 27.- Variation of relative reduction in induction time with 
equivalence ratio for = 1200 K and p = 4.0 atm. 
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Figure 28. - Variation of relative reduction in induction time with 
equivalence ratio for = 1200 K and p - 1.0 atm. 
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Figure 30.- Relative reduction in induction time as a function of effective mole fraction computed by 

empirical relationships for p = 1.0 atm and (j> = 1.0. 
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Figure 31.- Comparison of induction time for p = 1.0 atm and (^ = 1.0 using the kj 2 

and kj^2A ^ 10“^. 
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Reaction time, Ty, ^sec 



Figure 32.- Reaction time as a function of the reciprocal of mean temperature along 
isobars for 0 = 1.0 and Xq = 0 and 10-^. 
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Reaction time, r_, fxsec 
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Figure 33.- Reaction time as a function of initial temperature along isobars 
for (j) = 1,0 and Xq=s 0 and 10"^. 
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